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Abstract 

Using the Poisson bracket method, we construct the hydrodynamics of nu- 
clear matter in the chiral limit, which describes the dynamics of all low-energy 
degrees of freedom, including the fluid-dynamical and pionic ones. The hy- 
drodynamic equations contain, beside five Euler equations of relativistic fluid 
dynamics, Nf — 1 second order equations describing propagating pions and 
— 1 first order equations describing the advection of the vector isospin 
charges. We present hydrodynamic arguments showing that the pion velocity 
vanishes at the second order phase transition at Nf = 2. 
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Hydrodynamics is the theory describing the low-frequency, long-wavelength dynamics 
of liquids (or, in an extended sense, of any system). In this regime most degrees of freedom 
become irrelevant since they relax during the time characteristic of particle collisions; the 
only ones that survive are either those related to the conservation laws or the phases of 
the order parameters of broken continuous symmetries. The simplest example is normal 
fluids, where hydrodynamic variables arise due to the conservation of energy, momentum, 
and particle number. In superfluid He 4 and He 3 , additional hydrodynamic variables emerge 
from the symmetry breaking by the condensate. Although in all of these systems the physics 
is quite complicated at the molecular level, at large scales the hydrodynamic equations 
have simple forms dictated by the symmetries, the pattern of symmetry breaking, and the 
conservation laws. Such equations typically involve unknown coefficients, which can be 
computed from the microscopic theory or measured in experiments. 

A similar philosophy is shared by the chiral perturbation theory, which describes the 
long-distance dynamics of QCD with light quarks 0. At low energies, QCD is a strongly 
coupled theory where not much can be computed, at least at this moment, in a reliable 
fashion. However, well below the chiral scale (about 1 GeV), the dynamics is determined by 
the chiral Lagrangian, which can be written down knowing only the chiral symmetry and 
the pattern of chiral symmetry breaking of QCD. To the lowest order, pions are governed by 
the nonlinear sigma model, the only free parameter of which is the pion coupling constant 

which can be determined by matching the predictions of theory with experiment. 

In nuclear matter in the chiral limit, the low-energy degrees of freedom include both the 
fluid-dynamical variables (the energy-momentum tensor) and the chiral ones that describe 
the massless Goldstone modes arising from the breaking of chiral symmetry. All these degrees 
of freedom are coupled to each other; therefore, a full hydrodynamic treatment must include 
all these modes. In this respect, the hydrodynamics of nuclear matter is more similar to 
that of superffuids than of normal fluids. Treatments so far have largely dealt with the fluid 
dynamical and chiral variables separately, ignoring the interplay between the two || . 

The purpose of this paper is to construct the hydrodynamic theory of nuclear matter in 
the chiral limit, capable of describing all low-energy degrees of freedom of the latter. The 
primary place where such a theory can be applied is in the theory of heavy-ion collisions, 
where pion modes can arise by thermal activation or from the spinodal instabilities at the 
chiral phase transition. 

Degrees of freedom. — For definiteness, we will be working in Nf = 2 QCD, but most 
formulas remain valid at any value of Nf. As the first approximation, we take the quark 
masses to be zero, m u = rrid = 0. The nonzero pion mass, hopefully, can be taken care of by 
expanding around the chiral limit. We will assume that the temperature and the chemical 
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potential are low enough so that the chiral symmetry is broken. 

The first step toward a hydro dynamic description is to identify the hydrodynamic vari- 
ables. As mentioned above, the latter are the densities of conserved charges or the phases 
of order parameters of broken continuous symmetries. In nuclear matter, different vacua 
arising from chiral symmetry breaking are characterized by the phase £ of the condensate, 
where £ G SU(2), which transforms under chiral rotations as £ — > ISLE) . Slow variations 
of £ correspond to soft pions, which, due to the derivative coupling, relax much slower than 
a typical mode, and, hence, are hydrodynamic degrees of freedom. The conservation laws 
include those of the energy and the momentum, as well as the baryon number n = J dxq^q (q 
denotes quarks) , and the left- and right-handed isospin charges f<ixpL R (x) = /^x^l.rit^r, 
which are the generators of left and right isospin rotations. Therefore, the hydrodynamics 
variables are the energy and momentum densities T° v , the baryon number density n(x), the 
left and right isospin charge densities Pl iR (x), an d the SU(2) phase of the condensate £(x). 

Hydrodynamics of a reduced set of variables. — A nontrivial hydrodynamics can be already 
constructed for a reduced set of variables, chosen here to contain p£ R and £. In this 
simplified treatment, we ignore the variation of the fluid dynamical degrees of freedom, 
regarding the fluid as frozen. Such a theory is not the full theory, but its construction is 
much simpler and the equations obtained are suggestive of those in full hydrodynamics, 
which will be constructed later. 

The method we use to write down the hydrodynamic equations is the elegant Poisson 
bracket technique [§]||, which regards the ideal hydrodynamics of any fluid as a Hamiltonian 
system defined by a Hamiltonian and a set of Poisson brackets between hydrodynamic 
variables. The Poisson brackets arise from the canonical commutators in the microscopic 
theory; the ones between pl,r are the classical version of the current algebra commutators, 



while those between pl.r with £ are defined by the transformation laws of £ under chiral 
rotations, since pl,r are the densities of charges that generate these transformations, 



K(x), P l(y)} 
{pr(x), P R (y)} 

K(x), &(y)}, 







r 6c P L(x)5(x-y) 

r fec p R (x)5(x - y) 



(1) 



K(x),£(y)} = -zy£(x)5(x-y) 



K(x),£(y)} = i£(x) T 5(x-y). 



(2) 



Finally, the Poisson brackets of S with itself vanishes, 
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{E(x), E(y)} = 0. (3) 

Now let us turn to the Hamiltonian. We will limit ourselves to the leading order in 
derivatives of E, since we are interested in the dynamics at the largest scales. We will also 
keep only leading-order terms in powers of p. The most general form of the Hamiltonian 
consistent with chiral symmetry is 




where pl,r = Pl,r^~ an d / S) ft, and f v are constants with the dimension of mass, whose 
physical meaning will become clear later. 

Taking the Poisson brackets with the Hamiltonian @, we obtain the equations of motion 
of the hydrodynamic variables, 



<%E = -72 0°l E - Ep R ), (5) 
Jt 

dtPL = -diJt (6) 

dtPR = -diJ*, (7) 



where 

J\ = \f s W\ Jf = \f s ^d^. 

Equations (0) and (|7]) reflects the conservation of left- and right-handed flavor charges, where 
J 4 L ' R plays the role of the chiral isospin currents. Note that f v does not appear in Eqs. ©-(0); 
the reason is that tr(p L + £pR,£'0 2 = ^(Pl+Pr+^PlSprS^) is a Casimir operator. Equations 
(H)-(0) completely determine the dynamics of p and E. These differential equations are of 
first order in time and describe the evolution of 3(iV f 2 — 1) = 9 variables. 

The relation to the nonlinear sigma model. — The hydrodynamic equations ©-(0) should 
be contrasted with the field equations of the nonlinear sigma model, which describes the 
dynamics of pions at zero temperature and chemical potential. The latter is composed of 
three equations of second order in time derivative, which can be rewritten as six first order 
equations describing the evolution of E and c^E. Therefore our hydrodynamics contains at 
least three extra degrees of freedom that are not presented in the nonlinear sigma model. 

To find the relation of the hydrodynamic equations with the field equations of the non- 
linear sigma model, one solves Eq. (^) with respect to p^ and pr and expresses them via 
<9 t E and a new variable a, 

PL = -^ 2 S9 t Et + i«, (8) 
p R = -^ 2 Eta t E + ist a E. (9) 
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In particular, a = pl + SprE^, therefore in the vacuum where S = 1, a is the density of 
isovector charge. The equation of motion for a reads 

a = ~[Zd t Z\a], (10) 

while the equation for E is now second order in time derivative, 

if?d t {zd t ^) - ifidipatf) + a] = o. (ii) 

Equation (|T0|) always allows a = as a solution. In this case, Eq. fllTD reduces to the field 
equation of the nonlinear sigma model with pion velocity equal v n = f s j f t . Therefore / s and 
ft play the role of spatial and temporal pion decay constants, respectively || . The nonlinear 
sigma model can be interpreted as the Hamiltonian system ([]])-(§) with the constraint 
p L + £p R £t = o. 

In general, however, a needs not vanish. Equation (|H]) implies that, in the dissipationless 
limit we are considering, a only precesses with time. Since Eq. flH]| ) is first order in time 
derivative, once dissipation is included into the theory a will become a true diffusive mode. 
In any case, when a^O, our hydrodynamic equations are different from the field equations 
of the nonlinear sigma model. 

Hydrodynamic pion condensation. — To understand the physical meaning of the term 
proportional to a in Eq. (pTTJ) , consider the case where S makes only small variations around 
E = 1. As noted above, a is now the isovector charge. Subsequently, a can be nonzero if 
baryons are included into the theory. For example, if the baryon background contains more 
neutrons than protons, then a 3 = p p — p n is nonzero and negative. It is clear from the 
discussion above that a is the only baryonic degree of freedom that enter the hydrodynamic 
theory. 

By expanding around S = 1, using S = exp(z/^ 1 A a 7r a ), one finds the following linearized 
equation for the pion field on an isospin-asymmetric background: 

«9 V - vld^ + f t - 2 f abc d 7r b a c = 0, (12) 

Such an equation can also be obtained from the mean-field approximation of the chiral 
perturbation theory, by replacing in the interaction Lagrangian —^f 7 ^ 2 f abc 7i a d^7i b N'y fl Y^ 
the isospin baryon density Nj°yN by its mean value a c . We have shown that this mean- 
field procedure becomes reliable in the hydrodynamic limit, provided the in-medium pion 



decay constants ft and f s are used. Equation fll2| ) predicts a split between the dispersion 
relations of tc + and 7r~ in neutron-rich backgrounds where a 3 < 0. This is the hydrodynamic 
equivalence of pion condensation 0. 
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Pion velocity near the second order phase transition. — Let us consider N{ = 2 and assume 
that the baryon chemical potential p is zero or small enough so that the phase transition in 
temperature is second order ||. We will argue here that the pion velocity vanishes at the 
critical temperature. Indeed, as one approaches the critical temperature T c , the dependence 
of the Hamiltonian (§) on the phase of the condensate £ should become weaker, and, at 
T = T c , H should not depend explicitly on E. The latter can happen only if at the critical 
temperature f s = and f t = f v . Now since f v enters the Hamiltonian (|j) as a 2 / f 2 , it is 
related to the response of nuclear matter in the isovector channel. Such a quantity has no 
reason to vanish at the phase transition. Therefore one can expect f v 7^ at T c and, hence, 
ft also does not vanish at the critical temperature, in contrast to f s . Therefore the pion 
velocity v n = f s / ft approaches zero when T — > T c . 

Full hydrodynamics. — Now let us turn to the discussion of the full hydrodynamic equa- 
tions, which contains not only chiral variables but also the fluid dynamical ones. These 
equations can also be derived from the general Poisson bracket technique similar to the 
one used above. The hydrodynamic variables include, beside £ and Pl,r, the baryon den- 
sity n(x), the entropy density s(x), and the momentum density T ok (x). The nonvanishing 
Poisson brackets, beside those written in Eqs. (0)-(0), are 

{T 0i (x), A(y)} = A(x)c^(x - y), A = s, n, p L , R , 
{T°*(x), E(y)} = -$E(x)*(x-y), 



{T 0l (x), T ofc (y)} 



*(x-y). 



dx % KJ ' dy k . 

In particular, now tr(pL + ^Pr^) 2 is no longer a Casimir of the Poisson algebra: it has 
a nonzero Poisson bracket with T ok . The Hamiltonian is chosen in the most general form 
consistent with symmetries and containing only second order of pl,r and derivatives of S, 

H = J d*T 00 (x) =Jdx(e + fijtrdiZd^ + |tr(p L - ^ Pr Y)) 2 + ^tr(p L + Hp R ^) 2 
-2 Cfc tr(p L S9 fc S t + tfd k Ep R )), (13) 



where e, /y, a, a v , and Ck are functions of s, n and T ok . It is, however, more convenient to 
work with the conjugate variables: the temperature T = dT 00 /ds, the chemical potential 
fi = dT 00 /dn, and the velocity v k = dT 00 /dT ok . In the frame where v k = 0, e is determined 
by the nuclear equation of state, while fj = 5ijf 2 /4:, a = 2f~ 2 , a v = 2f~ 2 are functions of 
the local temperature T and chemical potential p, and = 0. The equations of motion for 
our variables can be found by taking the Poisson brackets with the Hamiltonian (|T3"D. The 
condition that <9oT 00 = —diT° l allows one to determine the velocity dependence of e, fy, a, 
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a v , and Ck, which turns out to be equivalent to the condition of boost invariance. The latter 
allows the final equations to be written in a relativistically covariant form, although our 
formalism is Hamiltonian in nature. After quite laborious, but straightforward, calculations 
||, one finds that the full set of hydrodynamic equations consists of (i) a second order 
equation for S, 

idM ~ fsWu'Hd^ + / s 2 £^£t] + K£<9 M £t, a] = 0; (14) 

(ii) a first order equation describing the advection and precession of the isovector charge a, 

d^a) = - l -[u^d^\a}- (15) 

(iii) the continuity equation for the baryon charge, 

d^vTn) = 0; (16) 

and (iv) the conservation of energy momentum, 

d,T^ = 0, (17) 

where the energy-momentum tensor T^ v is a sum of a fluid dynamical part and a field (pion) 
part, 

f 2 

T" u = (p + p)u»u v - pg^ + J -^ti{d^d v T) + d"£0"E + ). 

Moreover, both the energy density p and the pressure p receive a contribution from the pion 
field E and the density of isovector charge a. In particular, 

V = Po + \{ft ~ f 2 s )u^tid^d v ^ + ^tr^S^St + / -2 tra 2 ) (lg) 

where p = p (T,fi) is the pressure at E = 1 and a = 0, and p is related to p by Legendre 
transformation, 

p = po + ~(K + 1)(/, 2 - ffluWtTd^drf + (K — l)^tr^E^Et + + i)f^xa 2 , 

(19) 

where K = Td/dT + pd/dp. As in relativistic fluid dynamics, it can be shown that the 
conservation of entropy is a consequence of energy- momentum conservation (|TTD, baryon 
number conservation ([16]), and thermodynamic relations. The conservation of left- and 
right-handed isospin currents, 
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^ = ~\[{ft - ft)uWtfd v Z + / S 2 E^E] + lu^dE, 

follows from Eqs. ([[4]) and fll5|). 

Equations (14 ) and (|15|) are the straightforward generalization of Eqs. (|10D and (PI) . 
At small temperatures and chemical potentials, /t ~ f s and a ~ 0, the equation of motion 
for E decouples from the fluid flow and coincides with the field equation of the nonlinear 
sigma model. When E = 1, the hydrodynamic equations reduce to those of relativistic 
fluid dynamics. Thus, our equations include both the relativistic fluid dynamics and the 
nonlinear sigma model as special cases. This makes our equations ideal for the study of the 
evolution of disoriented chiral condensates [K| and any other situation where soft pions are 
in multiply occupied states. For example, the question of coupling between sound waves 
and pions can be addressed using these equations. 

The treatment of this paper could be improved in several directions. First, one can 
include the effects of quark masses by introducing the mass term tr(ME) + h.c. into the 
Hamiltonian, where M depends on T and /i. Because of the relatively large pion mass, 
such a term could be quite important. This will modify Eq. (|TJJ) for E and give additional 
contributions to the pressure and the energy density in Eqs. fll8|) and (|T9|), but the main 
features of the theory remain unchanged ||. Other next-to- leading terms in the chiral 
expansion can be also added. Second, the dissipative effects can be introduced, although 
one will certainly encounter the well-known problems of ambiguity and instability of viscous 
relativistic fluid dynamics fll |. 



In conclusion, let us notice that the method developed in this paper should be applicable, 
with some modifications, to the color-flavor-locking phase of finite-density QCD with iVf = 3 



The only difference is the breaking of U(l) baryon symmetry. The hydrodynamics 
contains one additional hydrodynamic variable [the U(1)b phase], which gives rise to the 
superfluid baryon number current. 
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